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ABSTRACT 


The  radiation  effects  of  high-energy  protons  and  electrons  on 
germanium  and  silicon  have  been  studied.  The  influence  of  the  damage  on 
the  Hall  effect  and  minority  carrier  lifetime  rroduced  by  1.  5-,  1.  3-, 

5-,  and  30-MeV  electons  is  compared  with  that  produced  by  30-  and  50- 
MeV  protons.  The  experimental  results  are  directly  compared  with  tne 
theoretical  predictions  of  the  damage.  It  is  concluded  that  the  higher- 
energy  recoils  resulting  from  the  proton  irradiations  are  more  effective 
in  producing  damage  than  those  resulting  from  the  electron  irradiations, 
based  on  the  predictions  of  the  total  number  of  the  defects  that  different 
types  of  particles  are  expected  to  produce.  The  reason  for  this  may  be 
that  the  higher-energy  recoils  produce  displacements  which  are  separated 
farther  and  thereby  influence  the  measured  properties  in  a  different 
manner. 

Complexities  in  correlating  proton  and  electron  damage  which 
arise  from  the  influence  of  impurities  and  defect  motion  during  irradia¬ 
tion  are  discussed.  A  proton  simulation  based  on  an  admixture  of  elec¬ 
trons  and  reactor  neutrons  is  discussed  theoretically. 


ftwicu*  png«  bl«J-  therefor*  ..et  rilmed 


CONTENTS 


I.  INTRODUCTION  . . .  .  .  1 


II.  THEORETICAL  CALCULATION  OF  DISPLACEMENT  EFFECTS 
FROM  PROTONS.  ELEC  IRONS.  AND  NEUTRONS  IN  SILICON 

AND  GERMANIUM  .  3 

2,  1  Introduction  . .  3 

2.  2  Primary  Recoil  Ionization  Energy  Loss  . .  6 

2.2  Rutherford  Scattering  . . 

2.4  Nuclear  Elastic  Scattering . 8 

2.5  Proton  Damage  . 10 

2.6  Electron  Damage .  20 

2.  7  Neutron  Damage .  22 

2.8  Conclusions .  28 


HI.  EXPERIMENTAL  PROCEDURES .  3  3 

3.  1  Sample  Preparation  .  3  3 

3.  1.  1  Sample  Cutting  and  Identification .  3  3 

3.1.2  Preparation  of  n-type  Silicon . . .  3  3 

3.  I.  3  Preparation  of  p-type  Silicon .  38 

3.1,4  Preparation  of  n-type  Germanium .  36 

3.  1.5  Preparation  of  p-type  Germanium  . .  36 

3.2  Hall  Effect  .  36 

3.  3  Carrier  Lifetime .  39 

3.  4  Dosimetry  . 45 

3.  4.  1  Beam  Calibration .  45 

3.4.2  Beam  Monitoring .  45 

3,5  Experimental  Fac ilities  .  48 


IV.  EXPERIMENTAL  RESULTS . . .  50 

4.  1  Hall  Sample  Irradiations .  50 

4.2  Carrier  Lifetime  Irradiations  .  62 

4.  3  Annealing  Studies .  64 


V.  DISCUSSION  .  68 

5.  I  n-Type  Silicon . 68 

5.2  p-Type  Silicon . 72 

5.  3  Germanium . 7  3 


VI.  CONCLUSIONS 


75 


REFERENCES 


77 


v 


Illustr  itione 


1.  Fraction  of  primary  recoil  energy  lost  in  displacement -type 

collisions . . .  7 

2.  Ratio  of  measured  elastic  scattering  cross  section  to  the 
calculated  Rutherford  cross  section  as  a  function  of 

center -of -mass  scattering  angle  for  30-MeV  protons .  9 

3.  Ratio  of  measured  elastic  scattering  cross  section  to  the 
calculated  Rutherford  scattering  cross  section  vs  center - 

of -mass  scattering  angle  for  aluminum . 11 

4.  Ratio  of  measured  elastic  scattering  cross  section  to  the 
calculated  Rutherford  scattering  cross  section  vs  center - 

of -mass  scattering  angle  copper . 12 

5.  Ratio  of  measured  elastic  scattering  cross  section  to  the 
calculated  Rutherford  scattering  cross  section  vs  recoil 

energy  for  silicon . 13 

6.  Ratio  of  measured  elastic  scattering  cross  section  to  the 
calculated  Rutherford  scattering  cross  section  vs  recoil 

energy  for  germanium . 14 

7.  The  number  of  atoms  displaced  by  recoil  atoms  with  energy 

AT  about  T  vs  T  for  silicon . 16 

8.  Total  nvmber  of  displaced  atoms  produced  by  recoil  atoms 

with  energy  from  T^  to  T  vs  T/Tjj  for  silicon . 17 

9.  The  number  of  atoms  displaced  by  recoil  atoms  with  energy 

AT  about  T  vs  T  for  germanium  . . IP 

10.  Total  number  of  displaced  atoms  produced  by  recoil  atoms 

with  energy  from  T^  to  T  versus  T  for  germanium . 19 

l  1.  Ratio  of  scattering  cross  section  to  Rutherford  cross  section 

vg  recoil  energy  for  electrons  of  energy  ^  4  MeV . .  21 

12.  TRIGA  neutron  spectrum  11  in.  above  reactor . .  23 

1  3.  Angular  distributions  used  for  neutrons  scattered  by  silicon, 


normalized  to  1  barn,  vs  the  cosine*  of  the  scattering  angle  .....  24 


14.  Angular  distributions  used  for  neutrons  scattered  by 

germanium,  normalised  to  1  barn,  vs  the  cosine  of  the 
scattering  angle . 25 

15.  Total  neutron  elastic  scattering  cross  sections  for  silicon 

and  germanium . . . . 2b 

16.  Sample  cutting  and  identification . 34 


vi 


17.  Sample  mounting  . . 37 

18.  Irradiation  chamber  for  samples . 38 

1Q,  Irradiation  sample  chamber  for  lifetime  measurements  .......  40 

20.  Sample  chamber  with  sample  mount  unplugged  . . 4! 

21.  Light  pulse  from  General  Badio  SUobotac . 4  1 

22.  Lifetime  instrumentation  diagram . 44 

23.  Cross  section  of  thin  calorimeter . 46 

2.4.  Carrier  removal  in  10 -ohm -cm  n-type  silicon  irradiated  with 

30 -MeV  electrons . SI 

25.  Carrier  removal  in  10 -ohm -cm  p-type  silicon  irradiated 

with  30 -MeV  electrons . 

26.  Carrier  removal  in  10 -ohm -cm  n-type  silicon  irradiated 

with  30 -MeV  protons . SI 

27.  Carrier  removal  in  10 -ohm -cm  p-typc  silicon  irradiated 

with  30 -MeV  protons . 54 

28.  Carrier  removal  in  0.  5 -ohm -cm  n-type  silicon  irradiated 

with  30 -MeV  electrons . SS 

29.  Carrier  removal  in  0.  5 -ohm -cm  n-type  i-ilicon  irradiated 

with  30 -MeV  protons . 56 

30.  Carrier  removal  in  0.  5 -ohm -cm  p-type  silicon  irradiated 

with  10 -MeV  electrons . 57 

31.  Carrier  removal  in  0.  5-ohm -cm  p-type  silicon  irradiated 

with  30 -MeV  protons . 58 

32.  Carrier  removal  in  l0-ohm-cm  n-type  germanium  irradiated 

with  30 -MeV  protons  . 59 

33.  Carrier  removal  in  10 -ohm -cm  p-type  germanium  irradiated 

with  30 -MeV  protons.  . . 60 

34.  Rate  of  change  of  reciprocal  lifetime  of  10-ohm-cm  n-type 

germanium  for  various  one rgy  radiations . 6  3 

35.  Carrier  annealing  in  10-ohm-cm  n-type,  floating -zone, 

refined  silicon . 65 

36.  Carrier  annealing  in  10-ohm-cm  n-type,  quartz  crucible, 

pulled  silicon. . 65 

37.  Carrier  annealing  in  0.  5-ohm-cm  n-type,  quartz  crucible. 

pulled  silicon.  . . 66 

vii 


38.  Carrier  annealing  in  0.  5-ohm-cm  n-type  germanium . 66 

39.  Carrier  annealing  ir.  10-ohm-cm  p-iype  germanium  . 67 

Tables 

I.  Theoretical  Predictions  of  the  Structure  of  the  Primary  Recoil 

Energy  Soectra  from  Different  Particle  Irradiations  . 30 

II.  Average  Current  Density  and  Particle  Flux  for  a  1-jiV/sec 

Calorimeter  Signal  . . 47 

III.  Carrier  Removal  Rates  and  Initial  Carrier  Concentrations . 61 

IV.  Lifetime  Degradation  Rates . 62 

V.  Comparison  of  Defect  Introduction  Rates . 69 

VI.  Measured  Property  Change  per  Calculated  Unannealed 

Displaced  Atom . . 


SECTION  I 


INTRODUCTION 

Radiation  produces  permanent  damage  in  semiconuuv.tor  materials 
by  displacing  the  lattice  atoms,  which  in  turn  may  form  electrically  active 
defect  structures  in  the  lattice.  These  defects  could  change  the  important 
electrical  properties  of  the  material,  such  as  the  conductivity,  carrier 
mobility,  number  of  carriers,  and  the  minority  carrier  lifetime,  and  thus 
the  characteristics  of  a  device  built  of  these  materials  are  altered  as  the 
result  of  irradiation. 

A  spectrum  of  electron  energy  irradiations,  1.5,  2.3,  5,  and 
30  MeV,  of  semiconductor  materials  are  compared  with  those  of  30-MeV 
and  50- MeV  proton  irradiations  in  order  to  determine  whether  proton 
damage  can  be  simulated  by  electron  damage.  The  effect  of  proton  dam¬ 
age,  of  course,  is  important  in  space  applications,  and  it  is  of  interest 
whether  the  eued  of  this  damage  is  the  same  as  for  electron-produced 
damage,  primarily  because  electrons  a  *e  a  more  available  radiation 
source.  The  experimental  results  are  directly  jmpared  with  theoretical 
predictions  of  the  damage  produced  by  several  types  of  irradiations.  From 
a  theoretical  standpoint,  a  comparison  of  30-MeV  electrons  with  protons 
is  of  primary  interest  as  they  are  most  likely  to  be  commensurate.  Also, 
comparisons  of  the  experimental  results  are  made  with  theoretical  pre¬ 
dictions  for  electrons  of  1.  5-,  2.  3-,  and  5-MeV  energy  and  for  protons 
of  50- MeV  energy.  These  low-energy  electrons  provide  a  check  on  the 
lower  limit  of  electron  energies  capable  of  simulating  protons  and  the 
50- MeV  protons  provide  a  check  on  the  effect  of  the  larger  displacement 
clusters. 

The  materials  irradiated  were  chosen  to  cover  a  spectrum  of 
semiconductor  materials.  The  twelve  chosen  consisted  of  0.  5-ohm-crn  and 
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1 0 -ohm -cm  samples  of  boron-doped  p-type  silicon,  phosphor js-doped 
n-type  silicon,  gallium-doped  p-type  germanium,  and  arsenic -doped 
n-type  germanium.  The  silicon  samples  consisted  of  both  floating-zone 
refined  and  quartz-crucible  pulled  materials  to  check  the  effects  ot  oxygen 
cGn*-  nt,  whereas  all  of  the  germanium  samples  were  quartz-crucible 
pulled  „ 

The  phy.  ical  properties  chosen  for  measurement  were  those  that 

are  most  intimately  „  asociated  with  device  operation  and  also  sensitize  to 

radiation:  the  carrier  concentration  and  carrier  lifetime.  The  carrier 

concentration  is  presented  as  a  removal  rate,  or  the  average  number  of 

» 

carriers  per  cubic  centimeter  removed  per  incident  particle  per  square 
centimeter,  and  the  carrier  lifetime  is  presented  as  a  rate  of  change  of 
inverse  lifetime  with  incident  particle  flux.  Post-irradiation  annealing 
experiments  were  made  in  an  attempt  to  identify  by  inference  some  of  the 
properties  of  the  defect  responsible  for  the  damage. 

A  unitary  positive  correlation  between  irradiations  with  different 
incident  particles  requires  a  constant  ratio  of  damage  rates  across  a 
spectrum  of  materials  and  damage  types.  Any  properties  inferred  about 
the  defects  responsible  for  the  damage  should  also  match.  However,  a 
lesser  degree  of  correlation  with  different  irradiations  would  still  be  use¬ 
ful  in  predicting  damage  rates. 
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SECTION  II 


THEORETICAL  CORRELATION  OF  DISPLACEMENT  EFFECTS 
FROM  PROTONS.  ELECTRONS,  AND  NEUTRONS 
IN  SILICON  AND  GERMANIUM 

2.  1  INTRODUCTION 

A  complete  theory  of  displacement  radiation  effects  would  describe 
the  sequence  of  events  occurring  to  produce  lattice  defects  and  would  relate 
the  character  of  the  resulting  defects  to  change  in  macroscopic  physical 
properties  of  the  irradiated  material.  The  sequence  of  events  described 
would  include: 

1.  Collisions  between  the  incident  radiation  and  the  target  to 
produce  recoil  atoms. 

2.  Motion  of  recoil  atoms  through  the  lattice  to  produce  subsequent 
generations  of  secondary  displacements,  leading  eventually  to 
the  dissipation  of  the  kinetic  energy  of  the  recoil  atoms  and  to 
*he  establishment  of  thermal  equipartition  of  kinetic  energy. 

3.  Subsequent  thermal -activated  motion  of  defects  to  form  secondary 
defect  structures  which  are  stable  at  t!.e  sample  temperature. 

4.  The  relation  between  the  microscopic  structure  of  the  defects 
and  the  macroscopic  physical  properties,  which  are  measured. 

It  is  obvious  that  such  a  complete  theory  is  not  available,  and  its 
lack  might  cause  one  to  despair  of  establishing  a  sound  theoretical  relation 
between  primary  irradiation  aad  physical  property  changes.  The  most 
imposing  obstacle  appears  to  be  a  description  of  thermal  motion  and 
secondary  defect  formation,  particularly  as  influenced  by  impurities  and 
defects  already  present  in  the  crystal  before  irradiation.  It  is  well  known 
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that  impurities  which  have  no  significant  effect  on  the  important  elect' ’.cal 
properties  of  a  semiconductor  may  influence  strongly  its  irradiation  response. 
Interstitial  oxygen  in  silicon  is  a  well-known  example. 

Fortunately,  in  attempting  to  establish  not  the  absolute  magnitude 
of  the  damage  but  only  the  correlation  between  damage  produced  by  dif¬ 
ferent  irradiating  particles,  it  is  possible  to  be  less  demanding  of  the  theory. 
Specifically,  it  is  obvious  that  if  each  of  two  irradiations  produces  the 
same  concentration  and  energy  spectrum  of  primary  recoil  atoms,  the 
resulting  damage  must  also  be  the  same.  Even  this  requirement  may  be 
too  severe.  A  high-energy  recoil  atom  tends  to  lose  most  of  its  energy 

by  a  large  number  of  collisions,  in  each  of  which  the  fractional  energy  loss 

(1) 

is  small.  Moreover,  the  higher-energy  recoils  {>100  keV  for  silicon) 
tend  to  lose  much  of  their  energy  by  ionization  until  the  residual  energy  is 
significantly  lower.  Hence,  as  far  as  the  displacement  structure  iu  con¬ 
cerned,  there  is  not  much  difference  between  a  100-keV  and  a  200-keV 
silicon  recoil  atom. 

If  it  is  also  assumed  that  the  physical  properties  measured  are 
sensitive  to  change  in  lattice  structure  over  regions  small  compared  to  the 
range  of  a  high-energy  recoil  atom,  it  can  be  asserted  that  the  damage 
should  be  a  function  only  of  the  total  energy  expended  in  displacement  pro¬ 
duction  (recoil  energy  minus  ionization  loss).  At  a  cursory  glance,  the 
secondary  recoils  produced  along  a  high-energy  primary  recoil  track  have 
an  energy  not  strongly  dependent  on  the  primary  recoil  energy,  and  hence 
their  number  increases  with  increasing  primary  energy,  but  the  character 
of  the  damage  may  not  change. 

This  discussion  suggests  that  there  is  an  excellent  chance  of 
establishing  a  correlation  between  two  irradiations,  provided  two  conditions 
are  met: 

1.  The  number  and  energy  spectrum  of  lower -energy  recoils  must 
be  the  same,  and 

2.  The  amount  of  energy  expended  in  displacement  production  by 

high-energy  recoils  must  be  the  same.> 
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The  transition  point  between  low-  and  high-energy  recoils  is  not  known,  but 
intuitively  it  should  correspond  to  displacement  clusters  containing  at  least 
tens  of  displacement  recoils  (i.  e.  ,  primary  recoil  energy  is  >1  keV). 

The  purpose  of  the  following  discussion  is  to  examine  in  detail  the 
energy  distribution  of  primary  recoils  resulting  from  electron,  proton,  and 
neutron  collisions  and  to  evaluate  the  possibility  of  correlating  electron  and 
proton  irradiations.  The  conclusion  to  be  noted  is  that  with  an  appropriate 
scaling  factor  for  the  Rutherford  cross  section,  high-energy  electrons  and 
protons  satisfy  the  first  condition  above,  but  that  for  this  scaling  factor 
influence,  the  total  displacement  energy  imparted  by  protons  via  high-energy 
recoils  is  larger  than  for  electrons.  This  extra  portion  could  be  adequately 
reproduced  by  adding  a  reactor  neutron  irradiation  to  the  electron 
bombardment. 

Protons  and  electrons  both  interact  with  an  atom  at  low-energy 
transfers  primarily  through  Rutherford  scattering,  and  the  scattering 
cross  sections  have  the  same  dependence  on  energy  transfer  for  both  at  the 
low  end  of  the  primary  recoil  energy  spectra.  For  protons  at  higher  pri¬ 
mary  recoil  energies,  nuclear  elastic  scattering  becomes  important,  and 
the  scattering  cross  section  results  in  a  dip  and  then  a  peak  in  comparison 
with  the  shape  of  the  Rutherford  scattering.  Further  oscillations  in  the 
cross  section  are  generally  less  important  as  these  occur  at  high-enough 
primary  recoil  energies  that  most  of  the  energy  of  the  recoil  is  lost  through 
electronic  collisions  and  only  a  small  part  goes  into  displacing  atoms.  The 
nuclear  elastic  effect  may  be  important  at  proton  energies  as  low  as  10  MeV 
and  must  be  taken  into  account  whenever  proton  damage  rates  are  considered. 

The  low-energy  part  of  the  proton  recoil  spectrum  may  be  simulated 
by  electrons  by  choosing  the  electron  energy  such  that  the  lessening  of  the 
electron  scattering  crose  section  at  high  primary  recoil  energies  due  to 
relativistic  effects  approximates  the  first  dip  in  the  proton  cross  section. 

The  following  peak  in  the  proton  cross  section  may  be  simulated  by  reactor 
neutron*  which  give  little  contribution  at  low  primary  recoil  energies  but 
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produce  apeak  at  higher  energies  similar  to  that  of  the  protons.  The 
quantities  of  electrons  and  reactor  neutrons  can  be  mixed  to  give  the  same 
total  number  of  displacements  as  protons  and  with  a  primary  recoil  atom 
energy  spectrum  similar  in  shape  up  to  muc.i  higher  energies  than  with 
electrons  alone. 

Nuclear  inelastic  scattering  was  neglected  in  calculating  the  proton 
damage  rates  and  it  is  possible  that  the  displacement  production  is  under¬ 
estimated  by  up  to  50  percent  in  germanium  at  50  MeV.  The  effect  is  less 
effective  for  silicon  and  lower  incident  proton  energies  but  still  may  be  as 
much  as  25  percent  in  silicon  at  30  MeV. 


2.  2.  PRIMARY  RECOIL  IONIZATION  ENERGY  LOSS 


As  the  energy  of  a  primary  recoil  atom  from  a  scattered  particle 
is  increased,  a  smaller  proportion  of  its  energy  will  go  into  displacement 
collisions,  a  greater  proportion  of  its  energy  being  lost  through  ionizing 
collisions  with  orbital  electrons.  The  fraction  of  energy  lost  in  elastic  or 
atomic  collisions  has  been  computed  by  Lindhard,^  assuming  a  Thomas- 
Fermi  potential  for  atomic  scattering  and  a  velocity-dependent  electronic 


stopping  power.  The  results  of  Lindhard's  calculations  for  silicon  and 
germanium  are  shown  in  Fig.  1  as  the  fraction,  f,  of  the  primary  recoil 
energy  that  goes  into  displacement-type  collisions.  This  fraction  is  used 
in  the  calculation  for  the  number  of  displacements  from  a  primary  recoil 
atom  instead  of  the  more  usual  arbitrary  cut-off  energy  beyond  which  a 
recoil  atom  is  assumed  to  produce  no  displacements. 

2.3.  RUTHERFORD  SCATTERING 


The  mechanism  by  which  charged  particles  impart  energy  to  a  lattice 

atom  is,  neglecting  nuclear  forces  and  screening  by  orbital  electrons, 

m 

Rutherford  or  Coulomb  scattering  with  a  differential  cross  section  of 


figure  l.  Fraction  of  primary  recoil  energy  lost  ’ n  c'jsplacement-type  collisions 


where  T  is  the  recoil  atom's  kinetic  energy,  is  the  maximum  energy 
transfer,  and  b  is  the  classical  distance  of  closest  approach  in  a  head-on 
collision. 

At  recoil  energies  much  greater  than  the  displacement  energy,  T^, 
the  number  of  displaced  atoms  per  primary  recoil  is  linearly  related  to  the 
recoil  energy  and  may  be  taken  as  0.  5  T/T^.  Using  the  number  of  dis¬ 
placed  atoms  produced  as  the  physically  significant  quantity,  the  number  of 
displacements  produced  by  primary  recoils  between  energies  Tj  and  is 


showing  that  each  logarithmic  energy  interval  contributes  equally  to  the 
number  of  displacements.  This  suggests  a  comparison  between  protons 
and  electrons  of  the  ratio  of  the  true  to  Rutherford  cross  section  versus 
the  logarithm  of  the  primary  recoil  energy. 


2.  4.  NUCLLAR  ELASTIC  SCATTERING 

Protons  and  neutrons  both  undergo  nuclear  elastic  scattering  in  the 

force  field  of  the  target  nucleus.  This  scattering  is  fairly  well  understood 

(4-7) 

in  terms  of  the  optical  model  of  the  nucleus,  and  a  number  of 
experiments  ^  have  been  performed  to  measure  the  scattering  cross 

section.  In  fact,  extensive  measurements  have  been  made  of  the  scattering 
of  protons  and  neutrons  on  aluminum  and  copper,  whereas  few  measure¬ 
ments  have  been  made  on  either  silicon  or  germanium.  Since  the  experi¬ 
mental  measurements  and  the  predictions  of  the  optical  model  both  show  a 
slow  variation  of  the  cross  sections  as  a  function  of  atomic  number,  the 
ernes  sections  for  aluminum  a.od  copper  have  been  used  for  silicon  and 
germanium.  The  scattering  cross  sections  for  30-MeV  protons  on  aluminum, 
silicon,  and  phosphorous  and  on  copper  and  gallium  (taken  from  Ref.  8), 
are  shown  in  Fig.  2  to  illustrate  the  slow  variation  between  adjacent 
elements.  The  ratio  of  the  measured  to  Rutherford  cross  section  is  plotted 
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to  accentuate  the  nuclear  elastic  effect  and  to  eliminate  the  elemental 
dependence  of  the  Rutherford  cross  section,  which  all  of  the  cross  sections 
approach  at  low  energy  transfers. 


2.  5  PROTON  DAMAGE 


For  incident  protons,  the  cross  section,  N^,  f°r  producing  dis¬ 
placements  from  primary  recoils  between  energies  Tj  and  is 


fT 2  T 

%  L  2T, 


-3T-mt<T>  dFdT  • 

d  R 


where  d<r/d<r^  is  the  ratio  of  the  measured  to  Rutherford  cross  section,  and 
f  is  the  fraction  of  the  recoil  energy  lost  in  displacement-type  collisions. 
The  Rutherford  cross  section  for  nonrelativistic  protons  is 

d<r  ,2  T  "Z2eJ 
R  _  wb  m  _ _ 0 

dT  4  _2  =  .  _  _2  ' 

T  AEqT 


where  A  and  Z  are  the  atomic  weight  and  atomic  number,  respectively,  of 
the  target  atom,  E^  is  the  incident  proton  energy,  e^  is  the  electronic 
charge. 


Ze  f  1  ♦ 


■«"*)' 


The  ratio  of  the  measured  elastic  scattering  cross  section  to  the 
Rutherford  cross  section  for  aluminum  and  copper  is  plotted  in  Figs.  3  and 
4,  respectively,  for  a  number  of  incident  proton  energies.  {The  references 
from  which  these  are  taken  ere  noted  on  the  figures. )  Because  of  the  smooth 
variation  in  the  cross  section  with  incident  proton  energy,  an  interpolation 
was  performed  between  40  and  57  MeV  to  give  the  50-MeV  cross  section. 

The  croes  sections  used  at  3C  and  50  MeV  for  silicon  and  germanium, 
together  with  the  function  f ,  are  plotted  in  Figs.  5  and  6,  respectively. 
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L'lastw  gcaiurring  t  roas  s«-  on  to  the  i  akulated  Rutherford 
i  t j on  vs  rci  oi.l  enerey  for  silicon 


Figu 
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For  the  purposes  of  numerical  calculation,  the  average  value  of 
and  f  are  taken 
this  reduces  Eq.  (1}  to 


d<r/dcr_  and  f  are  taken  over  the  interval  AT.  about  T  .  On  integration. 
R  1  i  ° 


N,  (T.)  =  ~L- 
up  l  dtr 

R 


f 


2  4 

it  Z  en  T  ,  .  , 
0  ,  i+l/2 

- In - 


T,  2 A  E  T  "*  T. 

rp  i  0  d  i  - 1  /  2 

i 


where 


T.  ,  =  T.  .  +  AT.  and  T  -  [T2  ,  +  T2  ,  I 

1+1/2  i-l/2  i  i  1  i+l/2  i~l/2J 


1/2 


to  gi’'e  the  aforementioned  logarithmic  energy  intervals.  Five  intervals 
were  taken  per  decade,  resulting  in 


ln  ,Ti+l/2^Ti-l/2*  ~  tl'4'2 


The  cross  section  for  producing  displacements  from  recoil  atoms 


with  energy  in  the  interval  ATi  about  is  then  given  by 


2  4 
e  —  ! 

N  (T.)  =0.231  — — 
dp  i'  AE0Td  doR 


f. 


m  T. 
T.  i 

t 


(2) 


and  the  total  cross  section  for  displacements  produced  by  recoil  atoms 
with  energy  from  to  T.,  (T.),  is 


72  4  Ti 
sZ  e, 


NJ  (T.)  =  0.  231  — 
dpt'  i'  AE 


0  < 
,T  L  d 


d<r 

do-. 


0  d  x  'R 


f. 


T 

T.  i 

i 


(3) 


Equations  (2)  and  (3)  are  plotted  for  silicon  in  Figs.  7  and  8, 
respectively,  and  are  plotted  for  germanium  in  Figs.  9  and  10,  respectively. 


for  protons  of  30-  and  50-MeV  incident  energy.  These  figures  are  plotted 

2  4 

in  units  of  0.  231  ttZ  e^/AE^  for  direct  comparison  with  the  recoil  spectra 
of  other  types  of  scattering.  Also  plotted  in  Figs.  8  and  10  are  the  results 
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Figure  10  Total  number  of  d'spiaced  atoms  produced  by  recoil  atoms  with  energy  from 
Tj  to  T  versus  T  for  germanium 


expected  from  pure  Rutherford  scattering  using  both  the  computed  value  of 
f  and  f  =  1  for  the  case  where  all  of  the  recoil  energy  is  assumed  to  be 
dissipated  in  displacement  type  collisions. 


2.6.  ELECTRON  DAMAGE 

The  electron  energies  of  5  and  30  MeV  considered  here  are  in  the 

(14) 

relativistic  rarsge,  so  the  calculations  of  McKinley,  which  are  plotted 
in  Fig.  11  ,  were  used  for  the  ratio  of  the  scattering  cross  section  to  the 
Rutherford  cross  section,  do7dcr  .  Then,  for  the  differential  scattering 
cross  section, 

.  ,  .2T 

dcr  _  dff  vb  m 

dT  =  d<r  4  2 

R  r 


where  for  the  relativistic  range  for  electrons,  b  is  contracted  by  y  = 

2  *1/2 

(1  *  p  )  in  which  p  =  v/c,  v  is  the  electron  velocity  and  c  is  the  velocity 
of  light,  and 


T 

m 


=  2 


E0tE0+2m0c  > 

MAc2 

A 


in  which  Eg  is  the  electron  kinetic  energy,  m^  is  the  electron  mass,  and 
is  the  mass  of  the  target  atom.  The  differential  energy  cross  section 
then  reduces  to 


dcr 

dT 
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Analogous  to  the  proton  case. 
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Equations  (4)  and  (5)  are  also  plotted  in  Figs.  7  through  10  for  5- 

2  4  2 

and  30-MeV  electrons.  These  are  given  in  units  of  0.  462  irZ  e^/M^c  for 
direct  comparison  with  the  primary  recoil  atom  energy  spectra  of  protons. 


2.7.  NEUTRON  DAMAGE 

The  neutron  spectrum  used  for  damage  calculations  is  that  calculated 
by  Larsen^ for  a  position  in  the  shielding  water  11  in.  above  the  TR1GA 
reactor.  The  logarithmic  energy  intervals  used  are  the  same  as  those  used 
for  the  previous  sections.  Figure  12  shows  the  spectrum  in  terms  of  neutrons 
in  each  energy  interval  used.  The  constancy  of  the  spectrum  at  lower 
energies  is  a  reflection  of  the  1/E  spectrum  assumed  at  energies  less  than 
10  keV. 

At  neutron  energies  much  above  a  few  hundred  keV,  the  scattering 
is  markedly  anisotropic,  so  it  was  necessary  to  use  the  angular  distributions 
of  the  elastic  scattering  given  in  Ref.  16.  The  resonances  particular  to 
separate  elements  were  smoothed  in  the  calculation  by  taking  rather  coarse 
energy  increments.  As  was  done  with  protons  and  for  the  same  reasons, 
the  measurements  on  aluminum  and  copper  were  used.  These  angular  dis- 
tributions  were  normalized  to  1  barn  and  interpolations  were  made  when 
necessary.  The  normalized  angular  distributions  used  for  each  neutron 
energy  interval  are  shown  in  Figs.  13  and  14  for  silicon  and  germanium, 
respectively.  These  normalised  angular  distributions  were  then  weighted 
by  the  total  cross  section  from  Ref.  17.  The  results  are  shown  in  Fig.  IS. 

The  cross  section  for  producing  a  recoil  atom  with  an  energy  AT. 
around  T^  from  a  scattered  neutron  of  energy  Er  is 
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Jittering  cross  sections  for  silicon  and  germr.nium 
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where 


dcr/dT|E  T  -nd  dcr/d«|E  T 
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are  the  differential  energy  and  angular  cross  sections,  respectively , 

A 

evaluated  at  E  and  T.,  and  T jr  is  the  maximum  energy  that  may  l  ' 
n  1  n 

imparted  to  a  recoil  atom  with  a  neutron  of  energy  E  . 

The  c/ats  section  for  producing  displacements  per  watt-stcord  of 

reactor  power  due  to  recoils  of  energy  AT.  around  T.  from  scattered 

neutrons  of  incident  energy  AF  around  E  ,  N,  (E  ,  T  ),  is  just  the  cross 

67  n  n  dn  n  l  J 

section  for  producing  a  recoil  atom  with  an  energy  AT^  around  T^  times  the 

number  of  neutrons  incident  with  energy  AE  around  E  per  watt- second  of 

n  n 

reactor  power  times  the  number  of  displacements  per  recoil  atom,  i.e.  , 
r  .  ,  .  r  T.  T 
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where  4K£n)AE^  is  given  in  Fig.  12  and  the  factor  of  10  converts  the 
cross  sections  from  barns  to  square  centimeters. 

Then  the  cross  section  for  producing  displacements  per  watt-second 
of  reactor  power  due  to  recoils  of  energy  AT^  about  T.,N^  (T^)f  using 
AT^  -  0.465  T^  for  our  logarithmic  energy  intervals  and  Tjr  =  (4/A) 

[i  +  (1/a)]’2  E/n,  is 
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where  the  summation  extends  over  the  reactor  neutron  spectrum  from  the 
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maximum  energy  E^  =  8  MeV  down  to  =  (A/4)[l  +  (1/A)  ]T^,  the  least 

energetic  neutron  that  can  give  energy  T^  to  a  recoil. 
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Also,  the  total  cross  section  for  displacements  produced  by  recoil 

atoms  with  energy  from  T  to  T.,  W,  (T),  is 

d  x  dnt  1 


VTi 


N,  T.  =0.730  x  10  4,4  A.  i  +  _  £  f  T2  V  AE  , 
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(8) 


which  is  analogous  with  that  for  protons.  Equations  (7)  and  (8)  are  plotted 
in  Figs.  7  through  10  in  units  such  that  the  total  number  of  displacements 
due  to  neutrons  is  equal  to  the  difference  in  the  total  number  of  displace¬ 
ments  produced  by  30-MeV  protons  and  30-MeV  electrons,  when  the  proton 
and  electron  displacement  production  rates  are  matched  at  low  reco:' 
energies. 


2.  8.  CONC  LUSIONS 

A  complete  discussion  of  the  theoretical  results  without  the  experi¬ 
mental  comparison  is  inappropriate  here  except  for  a  few  basic  conclusions. 
From  an  inspection  of  Figs.  7  and  9,  the  general  features  of  the  primary 
recoil  atom  spectra  from  electron,  proton,  and  reactor  neutron  collisions 
are  quite  clear.  Electrons  produce  displacements  at  a  nearly  constant  rate 
as  the  primary  recoil  energy  or,  equivalently,  the  displacement  cluster  si^e 
is  increased  up  to  nearly  the  maximum  energy  recoil.  The  maximum 
energy  recoil  is  approximately  proportional  to  the  square  of  the  incident 
electron  energy,  instead  of  being  linear  as  with  protons,  because  the 
relativistic  mass  of  the  electron  increases  with  energy.  For  protons, 
nuclear  elastic  scattering  produces  oscillations  in  the  displacement  pro¬ 
duction  cross  section.  The  positions  of  these  oscillations  are  nearly  constant 
in  recoil  atom  energy  as  the  incident  proton  energy  is  increased,  with  new 
oscillations  being  introduced  as  the  incident  proton  energy  and  maximum 
atom  recoil  energy  are  increased.  The  enhancement  of  the  scattering  cross 
section  due  to  nuclear  elastic  scattering  as  the  incident  energy  is  increased 
is  primarily  due  to  the  growth  of  the  first  one  or  two  peaks  in  the  scattering 
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cross  section.  The  effect  of  subsequent  peaks  is  reduced  by  the  function, 
f,  which  is  the  fraction  of  recoil  energy  going  into  displacement-type 
collisions.  From  an  inspection  of  Figs.  3  and  4,  it  can  be  seen  that  nuclear 
elastic  scattering  is  important  to  quite  low  incident  proton  energies, 
especially  in  silicon. 

When  30-MeV  electron  irradiations  are  compared  with  proton 
irradiations,  the  deviation  in  displacement  production  occurs  only  for 
displacement  clusters  of  ~500  displaced  atoms  for  germanium  and  ''-1000 
displaced  atoms  for  silicon.  When  5-MeV  electrons  are  considered,  the 
electron  simulation  is  much  less  likely  to  be  successful  since  the  clusters 
scale  only  to  10  and  20  displaced  atoms  for  germanium  and  silicon, 
respectively. 

Table  I  summarizes  the  results  of  these  calculations  with  a  dis¬ 
placement  energy  T^  equal  to  25  eV  for  silicon  and  32  eV  for  germanium. 

The  second  column  presents  the  relative  magnitude  of  the  low-energy  portion 
of  the  Rutherford  scattering  cross  section,  i.  e.  ,  a  number  by  which  the 
fluence  of  particles  should  be  divided  to  achieve  the  same  number 
and  spectrum  of  low-energy  recoils.  The  third  column  presents  the  recoil 
energy,  in  units  of  T^,  above  which  the  total  displacement  energy  deviates 
by  more  than  20  percent  from  the  Rutherford  value  for  large  maximum 
energy  transfer,  but  with  corrections  for  ionization  energy  loss.  This  quantity 
is  a  measure  of  the  range  of  validity  of  the  low-energy  recoil  spectrum 
simulation.  The  last  column  presents  the  effective  cross  section  for  pro¬ 
ducing  displaced  atoms,  irrespective  of  the  size  of  the  cluster  with  which 
the  atoms  are  associated.  Again,  this  number  should  be  used  only  for 
relative  magnitudes  since  most  of  the  displaced  atoms  undergo  thermally 
activated  rearrangement  at  room  temperature. 

It  is  apparent  from  the  foregoing  calculations  that  the  best  simulator 
for  proton  irradiations  to  meet  the  two  conditions  discussed  in  Section  2.  1 
is  a  combination  of  high-energy  electron  and  neutron  irradiation.  Specifically, 
the  following  correlations  are  indicated: 


In  silicon: 


2  2 
One  30-MeV  proton/cm  =  15.5  14-MeV  electron/crr.  + 

2 

0.  6 2  TRIGA  neutrons/cm 

One  50-MeV  proton/cm^  =9.3  14-MeV  electron/cm^  + 

C.  47  TRIGA  neutrons/cm^ 

In  germanium: 

2  2 
One  30-MeV  proton/cm  =  15.  5  20-MeV  electrons/cm  + 

0.  73  TRIGA  neutrons/crn^ 

One  50-MeV  proton/cm  =  9.3  20-MeV  electrons /cmr  + 

0.  70  TRIGA  neutrons/cm^ 


where  this  electron  energy  is  chosen  to  simulate  the  first  dip  in  the  proton 
nuclear  elastic  scattering  cross  section. 

In  lieu  of  irradiations  with  electrons  and  neutrons,  it  may  be  hoped 
that  less  accurate  but  satisfactory  simulation  is  possible  with  electrons 
alone,  usL.g  a  scaling  factor  based  on  total  number  of  displacements 
irrespective  of  primary  recoil  energy.  In  essence,  this  scaling  replaces 
some  of  the  higher -energy  proton -induced  recoils  by  low-energy  recoils. 

The  electron  simulation  will  then  overestimate  or  underestimate  the  damage, 
depending  on  whether  the  physical  property  measured,  taking  into  account 
thermal  motion,  is  more  or  less  sensitive  to  the  same  number  of  displaced 
atoms  being  distributed  in  smaller  recoil  clusters.  The  damage  ratios  on 
this  basis  are  as  follows: 


30-MeV  protons 
30-MeV  electrons’ 

50-MeV  protons 
30-MeV  electrons’ 

5-MeV  electrons 
30-MeV  electrons’ 

2.3-MeV  electrons 
30-MeV  electrons* 

1.  5-MeV  electrons 
30-MeV  electrons  ’ 


Silicon 

27 

18 

0.  58 

0.  42 

0.24 


Germanium 

20 

14 

0.  46 

0.  29 

0.  11 


It  should  also  be  noted  that  the  difference  in  the  contribution  of 
higher -energy  primary  recoils  due  to  nuclear  elastic  scattering  between 
30-  and  50-MeV  protons  is  found  to  be  less  than  the  enhancement  of  higher- 
energy  recoils  by  nuclear  elastic  scattering  over  Rutherford  scattering 
at  30  MeV.  It  was  found  that  for  silicon,  the  contribution  from  nuclear 
elastic  scattering  with  30-MeV  protons  is  twice  the  increase  in  nuclear 
elastic  scattering  when  going  from  30-  to  50- MeV  protons.  It  is  somewhat 
different  for  germanium  in  that  for  30-MeV  protons,  nuclear  elastic 
scattering  decreases  the  damage  calculated,  assuming  pure  Rutherford 
scattering,  by  approximately  the  same  amount  that  changing  from  30-  to 
50- MeV  increases  it.  The  50-MeV  irradiations  thus  only  serve  to  check 
results  obtained  at  30  MeV,  and  both  energies  will  be  affected  significantly 
by  nuclear  elastic  scattering. 
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SECTION  III 


EXPERIMENTAL  PROCEDURES 

3.1  SAMPLE  PREPARATION 

3.1.1  Sample  Cutting  and  Identification 

The  irradiation  samples  were  cut  from  the  semiconductor  boule  in 
the  manner  shown  in  Fig.  16.  First,  2-mm-thick  slices  were  cut  from  the 
boule;  they  were  then  remounted  and  cut  into  bars  1  mm  thick  by  20  mm 
long.  Of  the  six  slices  taken,  alternate  slices  (HI,  H2,  and  H3)  were  used 
for  the  Hall  coefficient  studies  and  the  other  slices  ( T 1 ,  r  2,  and  T  3)  were 
used  for  carrier  lifetime  studies.  Each  sample  was  identified,  as  shown, 
as  to  its  origin  and  eventual  use.  A  seventh  slice  was  taken  after  the  r3 
slice  and  cut  into  bars  for  the  unir  adiated  samples  required.  The  A^  and 
A^  pieces  outside  the  A^  samples  were  cut  1/2  and  1/4  mm  thick  to  pro¬ 
vide  samples  for  the  2.3  and  1.5-Mev  electron  irradiations,  respectively 

The  sample  bars  were  lapped  successively  with  400-grit,  600-grit, 
2/0,  and  4/0  emery  papers  and  then  they  were  etched.  The  germanium 
samples  were  etched  in  Super-oxal  etching  solution  (1  part  hydrofluoric 
acid,  1  part  hydrogen  peroxide,  and  4  parts  water),  whereas  the  silicon 
samples  were  etched  in  a  "white"  etching  solution  (equal  parts  of  glacial 
acetic,  nitric,  and  hydrofluoric  acids).  After  the  samples  were  etched, 
they  were  rinsed  in  distilled  water  and  electronic -grade  alcohol.  The 
preparation  of  the  sample  for  irradiation  then  varied  according  to  the  type 
of  sample. 

3.1.2  Preparation  of  n-Type  Silicon 

Each  n-type  silicon  sample  was  placed  in  a  quartz  jig  for  the  gold¬ 
bonding  operation.  The  jig  held  small  dots  of  0.5  percent  antimony -doped 


Scribe  line  along  length  of  boule 


Material  Boule  Resit  -  Pulled  30  MeV 

Length  tivity  Electron 


gold  against  the  sample  at  the  voltage  and  Hall  lead  positions.  The  jig  was 
spring  loaded  to  allow  for  contractions  due  to  the  melting  of  the  gold  dots 
during  the  bonding.  The  jig  and  sample  were  then  placed  in  an  oven  heated 
to  450  C  for  about  5  min  in  a  flow  of  forming  gas  and  then  they  were  cooled 
slowly.  After  cooling,  the  sample  was  removed  from  the  jig  and  the  ends 
for  the  current  leads  were  first  roughened  with  600-grit  emery  paper  and 
then  nickel  plated.  The  electroplating  solution  consisted  of  105  parts  by 
weight  of  nickel  ammonium  sulfate,  15  parts  of  ammonium  chloride,  and 
15  parts  boric  acid.  The  samp'e  was  then  washed  in  alcohol  and  acetone 
and  the  plated  ends  were  tinned  with  indium  solder. 

A  1-mil  copper  wire  was  soldered  to  each  gold  button  to  provide 
the  sample  lead.  A  copper-constantan  thermocouple  was  soldered  to  1 
gold  button  and  used  as  a  Hall  lead  on  the  Hall  coefficient  samples.  The 
thermocouple  was  attached  at  a  current  lead  for  the  carrier  lifetime 
samples.  Soldering  to  the  gold  buttons  was  accomplished  by  first  tinning 
the  copper  wire,  leaving  a  small  ball  of  solder  on  the  end  The  solder  ball 
was  then  heated  in  contact  with  the  gold  button  with  a  small  jet  of  hot  heli¬ 
um  gas.  Usually  a  little  Ruby  fluid  flux  was  needed  to  fuse  the  solder  to 
the  gold.  Current  leads  were  then  soldered  to  the  ends  of  the  sample, 
which  wa^  then  washed  in  alcohol  to  complet**  th*  preparation  of  the 
samcle. 

3.1.3  Preparation  of  p-Type  Silicon 

The  procedure  for  preparing  the  p-type  silicon  samples  differed 
from  that  used  for  n-type  silicon  in  that  boron-doped  gold  dots  were  used 
and  the  current  ends  were  plated  with  a  rhodium  plating  solution  obtained 
from  Englehard  Industries,  Inc.  Alternatively,  point  contacts  wert  some¬ 
times  used  for  Hall  and  voltage  probes.  These  were  phosphor -bronze 
wires  with  the  ends  nointed  and  coated  with  aluminum  that  were  spring 
loaded  against  the  sample  Contact  was  made  by  fusing  the  wire  to  the 
sample  by  passing  the  current  from  a  spark  coil  through  the  contact. 
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3.1.4  Preparation  of  n-T ype  Germanium 

For  n-type  germanium  the  Haii  an;,  voltage  leads  were  made  by 
first  placing  four  small  pieces  of  b0/40  soft  solder  on  the  sample.  The 
solder  pieces  were  held  in  place  by  the  surface  tension  of  the  Ruby  fluid 
flux  used  to  facilitate  soldering.  The  sample  was  then  placed  in  an  oven 
which  consisted  of  a  gla  ss  tube  wound  with  a  heater.  A  How  of  an  inert 
gas  was  maintained  through  the  glass  tube,  and  the  sample  was  heated  to 
the  melting  point  of  the  solder  to  bond  the  solder  dots  to  the  sample  The 
ends  of  the  sample  were  tinned  directly  with  indium  solder  for  the  current 
leads.  l«ads  were  attached  in  the  same  m ■?  as  ‘.hot;  for  n-tvpe  silicon. 

3 .  i  .  5  Preparation  of  p-Type  Germanium 

The  preparation  of  p-type  germanium  differed  from  n-tvpe  only  in 
that  indium  solder  was  used  for  all  c  >  facts 

3  l  HALL  EFFECT 

Samples  irradiated  for  measurement  of  carrier  removal  rate  were 
mo'mted  as  shown  in  Fig.  ??  vita,  tihosphor-bronre  spring  .  lips  holding  the 
samples  to  an  aluminum  block.  The  samples  were  electrically  insulated 
from  the  block  and  clips  with  Mylar  adhesive  tape.  The  block  was  then 
placed  in  a  water-cooled  chamber,  as  shown  in  Fig.  !H,  so  that  the  HI 
sample  was  nearest  the  lid  and  the  J{2  sample  was  in  the  middle  The 
electrical  leads  of  the  H2  simple  were  the  only  hi  ought  out  of  the 

box.  The  Hall  coefficient  of  the  H\  sample  only  was  monitored  during 
irradiation,  the  other  two  samples  were  measured  in  the  laboratory  before 
and  after  irradiation  as  a  check 

A  thermally  conductive  path  from  the  aluminum  block  to  the  sample 
chamber  was  formed  by  using  tin-foil  sW:ms  between  the  block  and  the 
bottom  and  lid  of  the  sample  chamber.  The  chamber  was  also  filled  with 
helium  gas  during  irradiation  for  more  efficient  heat  transfer. 
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Hall  coefficients  were  measured  in  a  field  of  4  kgauss  with  a 
current  of  about  2  ma  passing  through  the  sample.  The  sample  Hall  probes 
were  connected  to  a  microvolt  ammeter  through  a  voltage  bucking  box  to 
cancel  out  any  potential  difference  at  zero  magnetic  field  due  to  misalign¬ 
ment  of  the  probes.  The  meter  drove  a  chart  recorder  and  thus  a  record 
was  made  of  the  Hall  voltage  at  opposite  polarities  of  the  magnetic  field. 
The  difference  in  the  two  Hall  voltages  was  used  to  calculate  the  carrier 
density  in  order  to  cancel  out  ui)  contact  effects.  The  current  passing 
through  the  sample  was  measured  by  monitoring  the  voltage  drop  across  a 
1-kfl  series  resistor. 

Hall  coefficients  were  measured  at  the  irradiation  site  before  the 
samples  were  irradiated  and  the  irradiation  was  interrupted  at  approp¬ 
riate  intervals  for  further  measurements  until  the  Hall  voltage  had  changed 
by  25  percent,  at  which  time  the  irradiation  was  terminated.  Irradiated 
samples  were  then  taken  to  the  laboratory  and  there  the  Hall  coefficient 
was  measured  at  room  temperature  after  successive  5-min  anneals  at 
325°,  345°,  367°,  388°,  and  400°  K. 

A  temperature  record  of  the  H2  samples  was  made  during  the 
irradiations  and  the  beam  intensity  was  controlled  to  limit  the  temperature 
rise  to  a  few  degrees  Kelvin. 

3.3  CARRIER  LIFETIME 

Carrier  lifetime  samples  were  mounted  in  triple  mounts  in  much 
the  same  manner  as  the  Hall  samples  were,  but  because  the  rate  of  life¬ 
time  damage  was  generally  much  higher,  sample  heating  was  much  less  of 
a  problem  and  thus  plug-in  mounts  were  developed,  as  shown  in  Figs  19 
and  20.  The  miniature  9-pin  connectors  permitted  quick  sample  changing 
without  the  delicate  operation  of  manipulating  the  sample  leads  themselves. 
As  with  the  Hall  measurements,  only  one  of  three  lifetime  samples  was 
monitored  during  the  irradiation;  the  other  two  served  as  a  check  and  were 
measured  in  the  laboratory  before  and  after  irradiation. 
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Chamber  Tube 
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Figure  20.  Sample  chamber  with  sample  mount  unplugged 


The  carrier  lifetime  was  measured  by  observing  the  transient 
conductivity  of  the  sample  in  a  4-point  measurement  with  carriers  injected 
by  a  light  pulse  from  a  General  Radio  Strobotac.  the  General  Atomic  linear 
accelerator,  or  the  General  Atomic  flash  X  ray.  Figure  21  shows  the  shape 
of  the  light  pulse  at  two  different  repetition  rates  near  the  end  of  the 
instrument's  range.  The  high-rate  pulse  was  used  for  the  measurements 
because  of  the  shorter  decay  time,  even  though  the  peak  light  output  is 
much  less.  Thus  it  was  necessary  to  place  the  light  source  very  close  to 
the  sample. 

Provision  was  made  for  measuring  the  transient  current  through 
the  sample,  but  this  was  found  to  be  unnecessary  because  of  the  high 
impedance  source  used.  The  transient  voltage  was  measured  with  a  dc 
current  of  about  2  mA  in  the  sample.  Contact  effects  were  minimized  by 
using  the  difference  between  the  voltage  signals  taken  with  opposite  polarity 
currents.  Pictures  were  taken  of  the  transient  voltage  signal  on  an 
oscilloscope  camera  with  the  instrumentation  presented  in  block  diagram 
in  Fig.  22.  The  gain  of  the  preamplifier  system,  consisting  of  Tektronix 
Type  G  plug-ins  in  a  127  power  supplv,  was  remotely  controlled  by  means 
of  Ledex  stepping  motors  on  the  gain  switches  of  the  Type  G  preamplifiers. 
The  cable  drivers  are  capable  of  driving  the  93-ohm  cables  connecting  the 
irradiation  and  control  rooms  with  1-V  peak-to-peak  signals.  The  ovcr-ali 
sensitivity  of  the  system  is  0.  5  mV/cm. 

Transient  voltage  signals  thus  taken  were  assumed  to  be  inversely 
proportional  to  the  conductivity  (constant-current  assumption).  The  sig¬ 
nals  plus  the  dc  voltage  meatared  across  the  sample  with  no  carrier 
injection  were  read  into  a  computer  program  which  calculated  the  carrier 
lifetime  as  a  function  of  injection  levels.  From  these  data,  the  lifetime 
degradation  as  a  function  of  radiation  flux  was  determined.  During  electron 
irradiations,  the  lifetime  measured  with  the  Strobotac  was  checked  against 
the  lifetime  measured  when  a  short  (0.  1-psec)  pulse  of  relativistic  electrons 
was  used  to  inject  carriers.  Taese  measurements  were  made  to  calibrate 


42 


Repetition  rate  700  ppm 

Sweep  5  peec/ cm 


Repetition  rate  25,OOQ  ppn 

Sweep  0.5  psec/ cm 


ure  2i.  Light  pulse  from  General  Radio  Stro 


the  Strobotac  for  use  in  off-site  irradiations  where  no  other  carrier- 
injection  source  was  available.  Because  of  the  length  of  tail  on  the 
Strobotac  light  pulse  and  the  uneven  carrier  injection  in  the  sample, 
measurements  with  the  Strobotac  were  limited  to  lifetimes  in  excess  of  a 
n,  icrosecond,  with  the  calibration  being  required  when  the  lifetime  was 
below  about  lOpsec. 

After  irradiation,  the  samples  were  put  through  the  same  annealing 
schedule  as  the  Hall  samples,  with  the  lifetime  being  measured  with  the 
flash  X-ray  facility. 

3.  4  DOSIMETRY 

3.  4.  1  Beam  Calibration 

Calibration  of  the  irradiation  beam  was  accomplished  by  use  of  a 
thin  calorimeter  consisting  of  a  small  copper  block  suspended  in  a  con¬ 
stant-temperature  environment.  Figure  23  shows  a  cross  section  of  this 
calorimeter  with  the  copper  block  suspended  in  an  evacuated  water-cooled 
chamber  by  fine  copper-constantan  thermocouple  wires.  The  entrance  and 
exit  windows  are  1-mil  Mylar  affixed  to  the  chamber  with  epoxy. 

The  sensitivity  of  the  calorimeter  for  a  pulse  of  radiation  delivered 
in  a  time  short  compared  to  the  thermal  relaxation  time  of  the  calorimeter 
block  can  be  calculated  in  terms  of  the  energy  deposition.  For  the  case 
where  the  copper  block  is  initially  at  room  temperature,  the  dose.  D  (Cu). 
in  rads  (Cu)  is  related  to  the  thermocouple  voltage  increment.  AV.  by 

D  (Cu)  =  0.  975  x  109  AV  . 

Using  as  the  minimum  detectable  sig..al  a  voltage  change  of  1  pV/scc  or 
a  dose  rate  of  0.  975  x  10^  rads  (Cu)/sec,  the  minimum  average  c  urrent 
density  and  particle  flux  for  the  irradiations  considered  were  calculated; 
the  results  are  presented  in  Table  II.  All  fluxes  encountered  at  the 
irradiation  facilities  were  at  least  five  times  greater  than  the  minimum 
detectable. 
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Table  II 


AVERAGE  CURRENT  DENSITY  AND  PARTICLE  FLUX 
FOR  A  1-pV/SEC  CALORIMETER  SIGNAL 


Incident 

Radiation 

Average 

Current  Density 
(A/  cm^) 

Partic’e  Fbix 

(particles  cm^  sec) 

50-MeV  protons 

1.23  x  10'9 

7.69  x  109 

-10 

9 

30-MeV  protons 

8.34  x  10 

5.20  x  10 

30-MeV  electrons 

7.33  x  10*9 

4. 56  x  10 

5-MeV  electrons 

7  75  x  10‘9 

4.34  x  1010 

2.3-MeV  electrons 

8.06  x  10'9 

5.04  x  1010 

1. 5-MeV  elections 

8.14  x  io‘9 

5.08  x  1010 
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3.  4.  2  Beam  Monitoring 

For  continuous  monitoring  of  the  beam  during  irradiation,  a 
secondary-emission  monitor  was  used  for  both  the  electron  and  proton 
irradiations.  It  consisted  of  a  0.  008- in.  titanium  foil  in  an  evacuated 
chamber  v'ith  0.  001 -in.  titanium  windows.  The  monitor  was  placed  between 
the  beam  collimator  and  the  sample.  The  average  replacement  current  to 
the  foil  was  related  to  the  calorimeter  signal  when  the  calorimeter  was 
placed  in  the  sample  position  in  calibration  runs  before  and  after  each 
day's  irradiations.  During  the  irradiations,  the  replacement  current  was 
monitored  on  a  chart  recorder  and  this  record  was  used,  together  with  the 
calorimeter  calibration,  to  calculate  the  incident  particle  flux. 

3.  5  EXPERIMENTAL  FACILITIES 

The  5-MeV  and  30-MeV  electron  irradiations  were  performed  ut 
the  General  Atomic  Linac.  which  is  an  L-band  traveling-wave  electron 
accelerator.  It  can  produce  pulses  of  electrons  of  3  to  45  MeV,  pulse 
widths  varying  between  0.  01  ard  4.  5  psec,  repetition  rates  as  high  as 
720  pps,  and  beam  currents  up  to  1  A. 

The  low-energy  electron  irradiations  (1.  5  and  2.  3  MeV)  were  per¬ 
formed  at  the  Dynamitron  accelerator  located  at  the  Convair  Division  ol 
General  Dynamics  Corporation  i i.  San  Diego.  It  is  a  cascade- rectifier  type 
of  accelerator  capabli  of  currents  up  to  1  rnA  dc. 

The  proton  linear  accelerator  at  the  University  of  Southern 
California  was  used  for  the  initial  31.  5-MeV  proton  irradiations.  This 
accelerator  runs  at  a  frequency  of  300  Me.  hat  pulse  widths  of  450  psec. 
and  pulse  repetition  rates  of  1  5  to  20  pps. 

The  final  proton  experiments  were  run  at  the  Oak  Ridge  National 
Laboratory,  where  the  synchronous  cyclotron  was  used.  This  is  a  dc 
machine  capable  of  currents  in  the  milliampere  range;  the  energies  used 
were  35  and  56  MeV. 
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Lifetime  annealing  experiments  were  performed  at  the  General 
Atomic  flash  X-ray  facility,  which  produces  120-pser  pulses  of  600-keV 
X  rays  with  a  total  dose  of  about  1  r  per  pulse. 
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SECTION  IV 


EXPERIMENTAL  RESULTS 


The  experimental  accuracy  attributed  to  both  the  Hall  coefficient 
and  the  carrier  lifetime  sample  irradiations  is  about  20  percent,  including 
both  the  measurements  of  the  parameters  in  question  ana  the  determination 
of  the  amount  of  irradiation  the  samples  have  received. 

4.  1  HALL  SAMPLE  IRRADIATIONS 

Considering  the  Hall  effect  data  first,  the  typical  irradiation 
response  of  10-ohm-cm  phosphorus-  and  boron-doped  materials ,  both 
pulled  and  floating  zone,  is  shown  in  Figs.  24,  25,  26,  and  27  for  a  30-MeV 
electron  irradiation  and  a  30-MeV  proton  irradiation  and  a  30-MeV  proton 
irradiation,  respectively.  The  response  for  corresponding  0.  5-ohm-cm 
material  is  shown  in  Figs.  28,  29,  30,  and  31.  Similar  data  for  germanium 
are  shown  in  Figs  32  and  33,  respectively,  for  10-ohm-cm  arsenic-  and 
gallium-doped  material  during  a  30-MeV  proton  irradiation.  The  data  are 
summarized  in  Table  III  for  the  1.  5-,  2.  3-,  5-,  and  30-MeV  electron 
irradiation  and  for  the  30-  and  50-MeV  proton  irradiations. 

The  conductivities  measured  prior  to  irradiation  are  subject  to 
inaccuracies  owing  to  surface  and  contact  effects.  The  most  meaningful 
parameter  for  sample  identification  is  then,  in  view  of  the  radiation  response 
measurements  performed,  the  initial  carrier  concentration.  This  is  taken 
as  zero  irradiation  intercept  of  the  carrier  removal  curves.  his  procedure 
will  give  an  indication  of  the  initial  bulk  carrier  concentration  with  the 
effects  of  contacts  and  surfaces  negated  to  the  same  degree  as  they  are  in 
the  removal  rates.  Conductivities  may  be  calculated,  if  desired,  using 
known  values  for  the  mobility  in  silicon  and  germanium.  These  initial 
carrier  concentrations  are  also  given  in  Table  III. 
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gure  26.  Carr  er  removal  in  10-otim-cm  n-type  silicon  irradiated  with  30-MeV  protons 
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Table  III 


CARRIER  REMOVAL  RATES  AND  INITIAL  CARRIER  CONCENTRATIONS  . 


Irradiation  Energy 

|  1, 5-MeV  |  2.  i-MeV 

5-MeV 

SO -Me  V 

Material 

1  electrons  ]  electrons 

electrons 

p  rut  on  s 

Carrier  Removal  Rates 
(in  unit?  of  Carrie;  c/cm^  per  unit  flax) 


(i:j  units  of  carriers/cm^) 


Si  !  Or, 

6.  9x  101*1' 

fc.-ix  10U 

6.3x 1014 

8.  Ox  1014 

7.2x  10H 

6.8x  1014 

Si  lOnZ 

3,6x  i0H 

3. lx  JO14 

3.7x 1014 

4.2x  1014 

1 4 

4.  5x  10 

3. Ox  1014 

,  .  1 ' 

,  15 

,  15 

15 

1  5 

1  s 

Si  tOp 

I.  6 <  I If 

1.  6x  10 

1. 6x  10 

I.6x  10 

1.6x  10 

1.6x  10 

.  !  » 

15 

.  15 

15 

14 

1  s 

Si  10pZ 

1.9V  10 

1. 9x  10 

1. 9*  10 

1.9x  10 

1. 9x  10 

1.9x  10 

Ge  1  On  Z 

3.3x  10U 

3.  9x  1014 

4. Ox  1015 

3.  9*  1014 

-  i  ,,.14 

3.  6X  10 

3.  Ox  1014 

Ge  lOpZ 

2.9k  I014 

--- 

2,7x  1014 

2.  3x  1014 

1.  5x  !014 

1.6x  !014 

Si  0.  5n 

— 

— 

— 

0.9x  10H 

l.2.x  1016 

— 

Si  0.  5nZ 

— 

... 

1.6x  1016 

2.  2x  1016 

— 

Si  0.  5p 

... 

-  - 

... 

4.6x  1016 

4.  5x  1016 

... 

Si  0.  5pZ 

M  • 

— 

— 

5.2x  1016 

5.  lx  10lb 

... 

Ge  0,  SnZ 

as  m  at 

— 

•  •• 

5.8x  1015 

5.7X  1015 

— 

Ga  K 5pZ 

a  a  a 

... 

— 

15 

3.9x  10 

--- 

— 
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4,  2  CARRIER  LIFETIME  IRRADIATIONS 


The  compilation  of  the  lifetime  experimental  results  is  given  in 
Table  IV,  For  the  lifetime  simple*,  the  specified  resistivity  tolerance 
of  the  boules  of  ±20%  is  as  good  as  our  dc  measurements,  taking  into  con¬ 
sideration  the  inaccuracies  mentioned  in  Section  4.  1.  The  response  of 
10-ohm-cm  n-type  germanium  to  various  electron  energy  irradiations  is 
given  in  Fig.  34.  This  response  is  typical  of  the  type  of  results  obtained 
for  all  samples.  The  points  for  the  30-MeV  proton  irradiation  are  also 
shown  on  an  expanded  scale.  The  procedure  in  these  experiments  was  to 
measure  the  lifetime  of  the  samples  at  the  flash  X  ray  before  and  after  the 
experiments  that  were  not  done  at  General  Atomic.  In  these  ca-»es  the 
Strobotac  was  needed  to  measure  lifetimes  during  irradiation  and  the  flash 
X  ray  served  as  a  check  on  the  Strobotac  calibration  performed  with  similar 
samples  at  the  linear  accelerator. 


Table  IV 

LIFETIME  DEGRADATION  RATES 
(in  units  of  I0~^  cm*^aec~^) 


Irradiation  Energy 

Material 

1.  5-Mev 
electrons 

2.  3-MeV 
electrons 

5-MeV 

electrons 

30-MeV 
electr  3ns 

30-MeV 

protons 

Si  lOn 

3 

3 

7 

8 

600 

Si  10nZ 

0.5 

1.3 

2 

3.8 

100 

Si  lOp 

0.30 

0.45 

I.  5 

5.5 

500 

Si  IOpZ 

0.30 

0.45 

1.5 

5.5 

500 

Ge  IQnZ 

0.06 

0.20 

C.  20 

1.0 

27 

Ge  IOpZ 

0.025 

0.35 

0.06 

0. 14 

18 

eciprocal  lifetime  of  10-ohrr-cm  n-type  germanium 
radiations 


4.  3  ANNEALING  STUDIES 


Annealing  experiments  we  -e  performed  after  irradiation  by  holding 

the  sample  at  325°,  345°,  36 7°,  388°,  and  400°K  successively  for  5  min 

and  intermittently  cooling  to  300°K,  at  which  temperature  the  measure- 

ments  were  made.  For  those  samples  that  did  show  annealing,  the  results 

for  the  Hall  samples  are  shown  in  Figs.  35  through  39.  The  fraction  of 

the  carriers  recovered  is  defined  as  — - where  n„  is  the  initial  carrier 

n0  -  nf  0 

concentration,  n^  is  the  number  after  the  final  irradiation,  and  n  is  the 
number  after  the  previous  anneal.  No  annealing  was  seen  for  10-ohm-cm 
or  0.  5-ohm-cm  pulled  or  floating- zone  p-type  silicon.  0.  5-ohm-cm  floating- 
zone  n-type  silicon,  or  10-ohm-cm  n-type  germanium.  The  annealing  was 
the  same  within  20  percent  for  all  energy  irradiations.  For  the  lifetime 
samples,  with  one  exception  no  annealing  was  observed  for  any  sample,  but 
trapping  effects  were  evident  for  both  p-type  silicon  samples  and  the  floating- 
zone  refined  10-ohm-cm  n-type  sample.  For  p-type  silicon,  a  30  to 
50  percent  room  temperature  anneal  was  found  between  the  irradiation  site 
and  the  measurement  made  at  the  flash  X  ray.  Thus,  for  all  types  of 
radiaiums  considered  here,  the  annealing  data  show  no  evidence  of  the 
creation  of  different  types  of  defects  by  different  irradiating  particles. 
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CARRIERS  RECOVERED 


*E"C£"T  or  CARRIERS  RECOVERED  PERCENT  Of  CARRIERS  RECOVERED 


Figure  37. 


Carrier  annealing  in  0.  5-ohm-cm  n~type,  quarts  crucible, 
pulled  silicon 


Figure  38.  Carrier  Annealing  in  0.  5 -ohm* cm  n-type  germanium 


PERCENT  OF  CARRIERS  RECOVEf 


SECTION  V 


DISCUSSION 

It  is  now  possible  to  compare  the  experimental  data  with  theorv. 
Since  neutron  irradiations  were  not  included  in  the  experimental  program, 
it  is  not  possible  to  perform  a  direct  check  on  the  best  simulation  predic¬ 
tion  presented  in  Section  l.  8.  Instead,  a  comparison  between  experiment 
and  a  correlation  with  total  displaced  atoms  can  be  examined.  It  is  re¬ 
cognized  that  this  correlation  may  be  inaccurate  since  the  higher-energy 
and  higher-mass  particles  produce  a  larger  fraction  of  the  displacements 
via  n  ore  energetic  recoil  atoms.  This  comparison  is  presented  in 
Table  V  as  ratios  between  the  rate  of  introduction  of  defects  relative  to 
the  introduction  rate  in  the  same  material  by  30-MeV  electrons. 

Another  more  instructive  comparison  i6  achieved  by  calculating  the 
property  change  per  unannealed  displaced  atom  calculated  as  in  Section  11. 
These  data  are  summarized  in  Table  VI. 

The  radiation  response  presented  in  Table  VI  is  discussed  in  more 
detail  below. 

5.1  n-TYPE  SILICON 

The  lifetime  degradation  data  for  n-type  silicon  follow  the  theory 
fairly  well  for  the  electron  irradiations  from  1.5  to  30-MeV,  whereas  the 
30*MeV  proton  data  show  a  relative  increase  in  defect  production  over  theory. 
The  data  for  both  pulled  and  floating- zone  material  can  be  expected  to  show 
more  scatter  than  similar  irradiations  with  other  materials  because  of  a 
6-nisec  trapping  time  in  the  floating-zone  material  that  makes  an  accurate 
measurement  of  the  low  injection  level  difficult  and  the  unexpectedly  low 
unirradiated  bulk  lifetime  of  the  pulled  material.  The  latter  was  less 
than  1  psec,  which  is  an  order  of  magnitude  less  than  is  to  be  expected  for 
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MEASURED  PROPERTY  CHANGE  PER  CALCULATED  UNANNEALED  DISPLACED  ATOM 
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See  Table  V  for  explanation  of  svrr.bols. 


this  material.  This  not  only  affects  the  accuracy  of  the  data  (e.  g.  ,  room 
temperature  annealing  between  the  irradiation  site  and  post-irradiation 
measurements  at  the  flntui  M  ray  could  go  unnoticed  owing  to  the  difficulty 
of  measuring  such  a  short  lifetime  with  the  light  source')  but  also  tin  ap¬ 
plicability  to  other  materials  of  this  type,  sinc  e  a  ditierent  detect  than  is 
usual  may  be  producing  the  irradiation  response. 

These  data  may  be  compared  with  those  from  damage  to  p  on  n 

(18)  ( 1  9,  <J0) 

silicon  solar  cells  due  to  protons  and  dm*  to  electrons.  ’  The 
solar-cell  damage  measured  is  the  minority  carrier  diffusion  length, 
which  is  comparable  to  the  minority  carrier  lifetime.  The  electron  ir¬ 
radiations,  which  covered  bombarding  energies  from  0.4  to  40  Mo V, 
showed  substantial  agreement  with  simple  theory  over  the*  range  from  1.  s 
to  30  MeV.  However,  the  proton  damage  measured  with  17-  to  liQ-MeV 
particles  showed  a  slowly  falling  damage  rate  from  17  to  SO  MeV  that 
cannot  be  totally  ascribed  to  rising  contributions  from  either  nuclear 
elastic  scattering  or  star  production.  Above  SO  MeV,  the  damage  rate 
assumes  the  1  /E  character  expected  from  Rutherford  scattering,  Die 
floating-zone  data  are  in  agreement  with  this,  showing  thai  recoils  from 
30-MeV  protons  have  twice  the  predicted  effectivene f s  in  producing  damage 
than  recoils  from  l.S  to  30-MeV  electrons.  Also,  lifetime  degradation 
rates  measured  with  neutrons  from  a  TRIGA  reactor^  ^  give  a  property- 
change  per  unannealed  displaced  atom  of  8,  using  half  the  calculated  neutron 
cross  section  because  the  spectrum  used  in  the  irradiations  was  softer  than 
the  one  calculated  in  Section  U.  When  the  above  data  are  <  ompared  with 
the  theory  presented  in  Section  II,  they  are  consistent  with  an  assumption 
that  primary  recoils  of  the  order  of  10  ^  eV  are  approximately  four  times 
more  efficient  than  predicted  in  producing  damage  in  a  room-temperature 
irradiation  than  recoils  an  order  of  magnitude  lower  ir.  energy. 

For  the  floating-zone  material,  the  radiation  response  of  the  Hall 
coefficient  seems  to  follow  that  of  the  iitetimt  up  to  the  50-MeV  proton 
irradiation,  where  the  increase  in  damage  over  that  expected  levels  off 
instead  of  further  increasing  as  do  the  data  of  Rosenzweig  on  lifetime. 

n 


The  pulled  material  exhibits  a  gradual  increase  over  theoretical  damage 
as  the  irradiating  electron  energy  is  raised,  but  otherwise  follows  the 

behavior  of  the  floating -zone  material.  Carrier -removal  experiments 

(22)  (211 
with  neutron  irradiations  have  been  done  by  Stein'  'as  well  as  by  Curt  . s'  ' 

but  the*r  results  are  in  disagreement  by  factors  of  2  to  4,  with  Stein's 

removal  rates  bein^  the  higher.  Curtis's  results  give  a  property  charge 

per  unannealed  displaced  atom  of  0.  03  and  0.  025  for  floating-zone  and 

palled  silicon,  respectively,  which  are  much  lower  than  expected  in  the 

light  of  the  proton  nata  and  the  apparent  agreement  with  our  lifetime  data. 

5.2  p-TYFF  SILICON 

No  signific  nt  differences  between  pulled  and  floating-zone  p-type 
silicon  were  observed  in  either  lifetime  or  Hall  damage  rates  in  any 
irradiation. 

Lifetime  damage  may  again  be  compared  with  damage  to  n  on  p 

silicon  solar  cells  in  Refs,  18,  19,  and  21.  The  electron  irradiations 

from  0.4  to  40  MeV  show  results  much  like  those  derived  here,  with  the 

damage  rising  with  increasing  bombarding  energy  much  faster  than  is 

calculated.  In  fact  Carter  and  Downing  fit  their  data  to  the  square  of  the 

displacement  cross  section.  Our  electron  irradiations  fit  such  a  curve 

quite  »,..l,  but  the  30-MeV  proton  irradiations  show  a  smaller  increase  as 

do  «e  neutron  data  of  Curtis,  with  a  property  change  per  unannealed  dis- 

-8  2  -t 

placed  atom  of  7  x  10  cm  sec  ~.  The  proton  data  of  Rosenzweig  from 
1.35  to  130  MeV  shows  similar  behavior.  Below  8  MeV,  the  damage  rate 
rises  as  1/E  in  accordance  with  Rutherford  scattering,  but  above  10  MeV, 
the  damage  rate  falls  only  slightly  up  to  50  MeV,  where  it  resumes  an 
approximate  1/  E  behavior  again.  As  with  n-type  silicon,  calculations 
cannot  account  for  this  behavior,  even  when  all  the  recoil  energy  is  as¬ 
sumed  to  go  into  producing  displacements,  which  is  an  unlikely  condition. 
With  a  more  realistic  division  of  energy  between  displacements  and  ion- 

5 

Nation,  such  as  imposing  a  limit  of  about  10  eV  above  which  recoils  are 
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assumed  to  lose  all  their  energy  by  ionization,  the  divergence  from  the 
experimental  curve  is  much  greater  and  probably  cannot  be  ascribed  to 
star  production. 

The  Hall  data  follow  the  same  characteristics  as  the  lifetime  data, 

but  there  is  a  slower  increase  in  damage  efficiency  as  the  energy  rise-. 

This  behavior  is  also  indicated  in  data  on  carrier  removal  in  15-  to 

(19) 

1 50-ohm-cm  p-type  silicon.  The  30-  and  50-MeV  protons  seem  to 
follow  theoretical  predictions  in  producing  Hall  damage  in  p-type  silicon. 
Curtis’s  results  are,  as  before  with  n-type  silicon,  lower  than  expected, 
with  a  property  change  per  unannealed  displaced  atom  of  0.03. 

Thus,  in  both  n-  and  p-type  silicon,  the  higher-energy  primary 
recoil  atoms  seem  to  have  an  effect  that  scales  faster  than  the  total 
number  of  displaced  atoms  they  produce.  The  magnitude  depends  on  the 
material  and  indeed  even  or.  the  type  of  damage. 

5.3  GERMANIUM 

Lifetime  damage  in  n-  and  p-type  germanium  is  in  an  essentially 
constant  ratio  to  that  calculated  for  all  irradiations,  except  for  a  small 
rise  in  n-type  germanium  and  a  large  rise  in  p-type  germanium  with 
30-MeV  protons. 

Hall  irradiations  on  n-type  germanium  are  in  agreement  with  theory 
for  electrons,  but  again  protons  produce  damage  more  efficiently.  P-type 
germanium  is  a  special  case  as  the  removal  rates  change  with  carrier 
concentration  and  previous  irradiation  history,  which  makes  a  compar¬ 
ison  between  irradiations  hazardous. 

Again,  with  germanium  as  with  silicon,  it  seems  that  the  high- 
energy  recoils  have  an  effect  on  the  damage  that  rises  more  rapidly  than 
the  number  of  displaced  atoms. 

The  lack  of  differences  in  the  annealing  behavior  only  indicates 
that  the  particular  defects  which  anneal  in  the  temperature  range  300 *K  to 
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400  *K  are  formed  proportionately  by  all  of  the  irradiating  particles.  This 
observation  makes  no  inference  on  the  correlation  between  those  other 
defects  which  did  not  anneal  in  this  temperature  range. 


SECTION  VI 


CONCLUSIONS 


From  the  experiments  performed,  it  was  found  that  the  defect 


introduction  rate  rises  faster  with  rising  irradiating  energies  than  the 
number  of  unannealed  displaced  atoms.  It  is  apparent  that  the  effect  of 
displaced  atoms  on  both  lifetime  and  carrier  concentration  is  greater  if 
the  atoms  are  produced  by  higher-energy  recoils.  This  effect  is  magni¬ 
fied  further  if  it  is  noted  that  more  than  half  of  the  displacements  pro¬ 
duced  by  protons  are  associated  with  recoil  atoms  having  the  same 
energy  spectrum  as  those  produced  by  30-MeV  electrons.  Thus,  the 
apparent  increase  in  over-all  effectiveness  of  the  protons  must  be  attri¬ 
buted  to  less  than  half  of  the  damage.  Therefore,  an  apparent  increase 
of  damage  per  displaced  atom  of  a  factor  of  2  to  3  must  represent  an 
increased  effectiveness  of  displaced  atoms  produced  by  the  higher-energy 
recoil  atoms  of  a  factor  of  3  to  5.  Hence,  it  is  predicted  that  the  property 
change  per  unannealed  displaced  atoms  should  be  even  greater  for  neu¬ 
trons  than  for  the  highest  energy  protons  used,  and  indeed  this  is  the 


case  for  lifetime 


ers. 


(21,  22) 


damage  in  similar  materials  irradiated  by  other  work 


For  Hall  damage,  the  situation  is  less  clear  as  there  are  two 

(221 

pieces  of  conflicting  data,  but  if  the  higher  rates  of  Stein's  data'  'are 
correct,  then  the  above  conclusion  is  upheld. 

The  difficulty  in  interpreting  this  apparent  high  effectiveness  of 
primary  recoils  In  the  range  of  100  keV  in  producing  damage  is  height¬ 
ened  by  an  understanding  of  the  process  by  which  a  recoil  atom  loses  its 
energy.  ^  ^  The  majority  of  collisions  the  primary  recoil  suffers  in  slow¬ 
ing  to  an  energy  less  than  10  keV  (where  it  seems  not  to  produce  exces¬ 
sive  damage)  involve  small  energy  transfers  and  any  displacement 


clusters  produced  are  quite  far  apart  and  are  in  no  essential  way  differ¬ 
ent  from  those  produced  by  a  less  energetic  recoil.  However,  calcula¬ 
tions  show  that  a  1-kcV  recoil  produces  displaced  atoms  in  its  track 
about  4  lattice  constants  apart,  which  is  within  the  interaction  distance 
for  the  displacements  (4  to  6  lattice  constants),  whereas  a  100-keV 
recoil  along  the  initial  portion  of  its  track  produces  displacements  about 
15  lattice  constants  apart  along  with  scbstantially  more  ionization.  This 
combination  of  distantly  spaced  displacements  in  an  ionized  region  may 
well  produce  secondary  defects  in  different  proportions  from  the  dis¬ 
placed  atoms  from  lower  energy  recoils  which  find  themselves  fairly 
closely  spaced  in  a  neutral  region. 

It  would  seem  that  this  behavior  would  be  different  for  different 
materials,  and  for  different  defects  controlling  the  irradiation  response. 
This  does  seem  to  be  the  case  as  the  varying  results  between  material 
and  parameter  measured  show.  Also,  impurities  or  defects  not  involved 
in  the  defect  directly  responsible  for  the  radiation  response  may  affect 
or  even  control  the  difference  in  the  high-energy  recoil  damage  efficiency. 

Thus,  while  it  is  possible  that  the  above  data  can  be  used  for  pur¬ 
poses  of  proton  simulation  with  electrons,  extreme  care  must  be  taken 
to  ensure  that  the  defect  controlling  the  radiation  response  is  the  same 
as  that  measured  here  and  that  no  impurity  is  present  that  would  affect 
the  damage  rate  of  the  high-energy  recoils  present  with  protons.  A  sim¬ 
ulation  with  a  better  chance  of  success  would  be  an  admixture  of  electrons 
and  reactor  neutrons,  as  described  in  the  theoretical  portion  of  this  report. 
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IS  AnTAACT 


The  radiation  effects  of  high-energy  protons  and  electrons  on 
germanium  and  silicon  have  been  studied.  The  influence  of  the  damage  on  the 
Hall  effect  and  minority  carrier  lifetime  produced  by  1.  5-,  2.  3  ,  5-,  and  30-MeV 
electrons  is  compared  with  that  produced  by  30-  and  50-MeV  protons.  The 
experimental  results  are  directly  compared  with  the  theoretical  predictions  of  the 
damage.  It  is  concluded  that  the  higher-energy  recoils  resulting  from  the  proton 
irradiations  are  more  effective  in  producing  damage  than  those  resulting  from  the 
electron  irradiations,  based  on  the  predictions  of  the  total  number  of  the  defects 
that  different  types  of  particles  are  expected  to  produce.  The  reason  for  this  may 
be  that  the  higher-energy  recoils  produce  displacements  which  are  separated 
farther  and  thereby  influence  the  measured  properties  in  a  different  manner. 

Complexities  in  correlating  proton  and  electron  damage  which  arise 
from  the  influence  of  impurities  and  defect  motion  during  irradiation  are  discussei 
A  proton  simulation  based  on  an  admixture  of  electrons  and  reactor  neutrons  is 
discussed  theoretically. 
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